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The driven lattice gas (DLG) evolving at low temperature helps understanding the kinetics of 
pattern formation in unstable mixtures under anisotropic conditions. We here develop a simple 
theoretical description of kinetics in Monte Carlo simulations of the DLG. A Langevin continuum 
analog is also studied which is shown to exhibit the same behavior. We demonstrate that pattern 
growth is mainly a consequence of single-particle processes and that, after a short transient time, 
in which a surface evaporation/condensation mechanism is important, hole diffusion in the bulk 
becomes dominant. Consequently, there is a unique relevant length that behaves I it) ~ t 1 ^ 3 for 
macroscopic systems except at some very early (perhaps unobservable) time. This implies sort of 
self-similarity, namely, the spatial pattern looks alike, but for a (non-trivial) change of scale at 
different times. We also characterize the structure factor, in which we identify Guinier and Porod 
regions, and its scaling form with both time and size. The underlying anisotropy turns out to be 
essential in determining the macroscopically-emergent peculiar behavior. 
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INTRODUCTION 

Many alloys such as Al-Zn, which are homogeneous 
at high temperature, undergo phase separation after a 
sudden quench into the miscibility gap (for details, see 
the reviews 0-0, for instance). One first observes nu- 
cleation in which small localized regions (grains) form. 
This is followed by "spinodal decomposition". That is, 
some grains grow at the expense of smaller ones, and 
eventually coarsen, while their composition evolves with 
time. In addition to theoretically challenging, the details 
are of great practical importance. For example, hardness 
and conductivities are determined by the spatial pattern 
finally resulting in the alloy, and this depends on how 
phase separation competes with the progress of solidifi- 
cation from the melt. 

A complete kinetic description of these highly non- 
linear processes is lacking. |5( Nevertheless, the essential 
physics for some special situations is now quite well un- 
derstood. This is the case when nothing prevents the 
system from reaching the equilibrium state, namely, co- 
existence of two thermodynamic phases. The simplest ex- 
ample of this is the (standard) lattice gas evolving from a 
fully disordered state to segregation into liquid (particle- 
rich phase) and gas (particle-poor phase) . (Alternatively, 
usingthe language of the isomorphic lattice binary al- 
loy, |fj the segregation is into, say Al-rich and Zn-rich 
phases.) As first demonstrated by means of computer 
simulations, 0, this segregation, as well as similar 
processes in actual mixtures exhibit time self-similarity. 



This property is better defined at sufficiently low temper- 
ature, when the thermal correlation length is small. The 
system then exhibits a single relevant length, the size £ (t) 
of typical grains growing algebraically with time. Con- 
sequently, any of the system properties (including the 
spatial pattern) look alike, except for a change of scale, 
at different times. 

This interesting property is revealed, for example, by 
the sphcricalizcd structure factor S (k, t) as observed in 
scattering experiments. After a relatively short tran- 
sient time, one observes that S (k, t) ~ J (t) • F[k I (t)} . 
Taking this as a hypothesis, one may interpret J and 
£ as phenomenological parameters to scale along the S 
and k axes, respectively. The hypothesis is then widely 
confirmed, and it follows that J (t) ~ £ (t) d where d 
is the system dimension. It also follows that F(x) = 
• \& (<7x) where 4> and 4" are universal functions. In 
fact, <E> describes the diffraction by a single grain, ^ is a 
grain interference function, and a characterizes the point 
in the (density— temperature) phase diagram where the 
sample is quenched. It then ensues that iff ~ 1 except at 
small values of k, so that, for large x, F(x) becomes al- 
most independent of density and temperature, and even 
the substance investigated. 0, 

The grain distribution may also be directly monitored. 
A detailed study of grains in both microscopy experi- 
ments and computer simulations confirms time scale in- 
variance. More specifically, one observes that the relevant 
length grows according to a simple power law, £ (t) ~ t a , 
and one typically measures a — 1/3 at late times. This 
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is understood as a consequence of diffusion of monomers 
that, in order to minimize surface tension, evaporate from 
small grains of high curvature and condensate onto larger 
ones {Ostwald ripening). In fact, Lifshitz and Slyozov, 
and Wagner independently predicted £ ~ t which is 

often observed, even outside the domain of validity of the 
involved approximations . 1 1 Ojj In some circumstances, one 
should expect other, non-dominant mechanisms inducin. 
corrections to the Lifshitz-Slyozov- Wagner one.Q, 0, 
For instance, effective diffusion of grains (Smoluchowski 
coagulation) leads to a — 1/6, which may occur at e arly 
times: [llT| interfacial conduction leads to a = l/4:|I3,ll3| 
and, depending on density and viscosity, a fluid capa- 
ble of hydrodynamic interactions may exhibit crossover 
with time to viscous (a = 1) and then inertial (a = 2/3) 
regimes. 4] 

Extending the above interesting picture to more re- 
alistic situations is an open question. The assumption 
that the system asymptotically tends to the coexistence 
of two thermodynamic (equilibrium) phases is often un- 
justified in nature. This is the case, for example, for 
mixtures under a shear flow, whose study has attracted 
considerable attention, e.g.Q|"E3- The problem is that 
sheared flows asymptotically evolve towards a nonequi- 
librium steady state and that this is highly anisotropic. 
Studying the consequences of anisotropy in the behav- 
ior of complex systems is in fact an important challenge 
(see, for instance, ^^-|13)- Another important exam- 
ple is that of binary granular mixtures under horizontal 
shaking. The periodic forcing causes in this case phase 
separation and highly anisotropic clustering. |42j 

In this paper, we study in detail the kinetics of the 
driven lattice gas (DLG) [13] following a deep quench. 
Our motivation is twofold. On one hand, the DLG is rec- 
ognized to be an excellent microscopic model for nonequi- 
librium anisotropic phenomena. [2l| On the other, the 
DLG is not affected by hydrodynamic interactions, which 
makes physics simpler. Our goal is timely given that the 
asymptotic state of the DLG is now rather well under- 
stood, and previous studies of kinetics altogether reveal 
an intriguing situation. HJ-H^ Following this pioneering 
effort, we here present a new theoretical description of 
the essential physics during anisotropic, noncquilibrium 
pattern growth. This is compared with new extensive 
computer simulations. A brief and preliminary account 
of some of our results was presented elsewhere. |2^ 



MODEL AND SIMULATION DETAILS 

The DLG consists of a d— dimensional, e.g., simple- 
cubic lattice with configurations n = {rij; i = 1, N}. 
The variable at each lattice site has two possible states, 
iii = 1 [particle) or [hole). As for the standard 
lattice gas, dynamics is a stochastic process at tem- 
perature T consisting of nearest-neighbor (NN) parti- 



cle/hole exchanges. This conserves the particle density, 
p = N^ 1 J2i n i: an d depends on n. 

A distinguishing feature of the DLG is that exchanges 
are favored in one of the principal lattice directions, say x. 
Therefore, assuming periodic (toroidal) boundary condi- 
tions, a net current of particles is expected to set in along 
x. This is accomplished in practice by defining a biased 
transition rate. We shall refer here to the energy function 
H = —4,JJ2nn n i n ji which describes attractive interac- 
tions between particles at NN sites, and to the transition 
rate (per unit time):[2l| 
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n* stands for configuration n after jumping of a particle 
to a NN hole; AH = H (n* ) — H (n) is the energy change 
brought about by the jump; and units are such that both 
the coupling strength J and the Boltzmann constant are 
set to unity. One further defines 8 = (^1,0) for NN 
jumps along ±x or along any of the transverse directions, 
say y, respectively. Consistent with this, E = Ex may 
be interpreted as a field driving particles, e.g., an electric 
field if one assumes that particles are charged. (One may 
adopt other interpretations, e.g., the binary alloy one. [a] 
Dynamics then consists of interchanges between particles 
of different species, one of them favored along x.) 

The DLG was described as modelling surface growth, 
fast ionic conduction and traffic flow, among a number 
of actual situations of practical interest. [2lJ A common 
feature in these situations is anisotropy, and that steady 
states are out of equilibrium. Both are essential features 
of the DLG induced by the rate Q). The only trivial 
case is for E = 0, which reduces to the Metropolis 
algorithm. In this case, detailed balance holds, and one 
simply has the familiar lattice gas with a unique (equi- 
librium) steady state. For any, even small E, qualita- 
tively new behavior emerges. In fact, detailed balance 
breaks down for E > and, consequently, the steady 
state depends on o;(n — > n*). Increasing E, one eventu- 
ally reaches saturation. That is, particles cannot jump 
backwards, i.e., —x, which formally corresponds to an 
infinite field {E — oo). 

The way in which the microscopic anisotropy Q con- 
veys into macroscopic behavior is amazing. 21] Con- 
sider, for simplicity, d — 2, p — | and E = oo. The 
system then exhibits a critical point at T = T^f ~ 
L4T C [E = 0), where TJE = 0) « 2.2691, with novel 
critical behavior. [13, Hg] Furthermore, the asymptotic, 
steady states below T^f do not comprise equilibrium 
phases. Instead, one observes a particle current and 
fully anisotropic phases; both are nonequilibrium fea- 
tures. The intensity of the current increases with T, 
and suddenly changes slope at T<2P (in fact, this prop- 
erty may serve to accurately locate the critical point). 
The stable ordered configurations consist of one stripe, to 
be interpreted as a liquid (rich-particle) phase of density 
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Pl {T) . The gas (poor-particle) phase of density pa (T) 
fills the remainder of the system. Except for some micro- 
scopic rou ghn ess, the interface is linear and rather flat, 
in general. [46j 

The computer evolutions reported here always begin 
with a completely disordered state to simulate the system 
at infinite temperature. We then model a sudden quench 
and the subsequent time evolution. With this aim, one 
proceeds with rate QJ that involves the temperature T 
at which the system is quenched. The run is followed 
until one stripe is obtained (eventually, in order to save 
computer time, the run was sometimes stopped before 
reaching the final stationary state). The code involves a 
list of -q it) particle-hole NN pairs from where the next 
move is drawn. Time is then increased by At = r\ (i) 1 , 
so that its unit or MC step involves a visit to all sites on 

the average. 

The lattice is rectangular, L\\ x L±, with sides ranging 
from 64 to 256 and, in a few cases, 512. Results concern 
an average over around thousand independent runs. Due 
to the great computational effort which is consequently 
involved, this paper describes simulations concerning a 
single point of the two-dimensional DLG phase diagram. 
That is, most of our evolutions are for p = A and E = oo, 
and simulate a quench at T = 0.8T C (E = 0) ~ 0.67^. 
This choice is motivated by the fact that clustering is 
then reasonably compact, which helps to obtain good 
statistics, while it proceeds fast enough, so that one can 
observe full relaxation to the steady state. In spite of this 
restriction, brief investigation of other points, together 
with some of our observations below, led us to believe 
that the validity of our results extends to a large domain 
around the center of the miscibility gap; in fact, such 
generality of behavior has been reported for E = 0.0,0, 

Sua 

GROWTH OF ORDER 

The DLG exhibits different time regimes during phase 
separation. Though they parallel the ones for E = 0, 
the peculiarities induced by the anisotropic condition are 
essential. 

Starting from complete disorder, there is a very short 
initial regime in which small grains form. The novelty 
is that typical grains are now fully anisotropic, stretched 
along x. The grains then rapidly coarsen to form macro- 
scopic strings, as illustrated in fig.l. Sheared fluids 
(an experimentally accessible situation that also involves 
both nonequilibrium physics and anisotropy) seem to ex- 
hibit similar initial regimes. 0, 0] That is, during a 
short time interval, they show larger growth rate along 
the flow than in the other directions, which is assumed 
to correspond to the initial formation of anisotropic re- 
gions. Afterwards, sheared fluids develop string-like 
macroscopic domains similar to the ones in the DLG. 



Fig.l includes a comparison with the zero-field case, 
i.e., the standard, isotropic lattice gas (LG). This clearly 
illustrates the strong anisotropy of nucleation and early 
phase separation for the DLG. Close inspection of these 
and similar graphs also seems to indicate relatively small 
but significant differences in the degree of segregation be- 
tween the two cases at a given time. That is, at small dis- 
tances, there is a more homogeneous distribution of par- 
ticles, both longitudinally and transversely, in the DLG 
than in the LG. The latter shows up more segregated at 
the same time, which is already rather evident by direct 
inspection of graphs for 1 < t < 100 in fig.l. We believe 
this reveals the different role played by surface tension as 
the degree of anisotropy is varied: Typical DLG grains 
are rather linear except at their longitudinal ends, where 
curvature may be even stronger than for the spherical 
clusters in the LG at comparable times. This seems to 
be at the origin of a smother transverse distribution of 
particles in the DLG at early times. On the other hand, 
the field also tends to smooth things longitudinally. 

In order to quantify the aforementioned observation, 
we evaluated the number of broken bonds in the direction 
of (perpendicular to) the field, n\\{t) (n±(t)) as a function 
of time during the early evolution stage. Then A(t) = 
[n±(t) + n.n (t)]/2JV is the density of broken bonds. The 
higher the degree of segregation at time t, the smaller is 
A(t). For instance, we observe in a large 256 x 256 lattice 
that A(t = 10) = 0.295 and A(t = 10) = 0.38 for the LG 
and DLG respectively, confirming the above observation. 
On the other hand, let B(t) = [n±(t) - n\\(t)]/2N. One 
would expect B(t) » (up to fluctuations) only for the 
isotropic system. In fact, we measured B(t) w for 
the LG, while B{t) rapidly converges to a nonzero value 
B(t) w 0.05 for the DLG at early times (again for a large 
256 x 256 lattice). We take this number, B(t) — 0.05, as 
characterizing the anisotropic shape of DLG clusters at 
early times. 

The difference of segregation between the DLG and the 
LG at early times merits further study. This will need 
to take into account the anisotropy of surface tension. In 
any case, this concerns a regime very near the initial, melt 
state that only bears minor practical importance, given 
that it extends extremely shortly on the macroscopic time 
scale. We are interested in the rest of this paper on the 
subsequent evolution, to be described on the assumption 
of a simple flat interface, which holds in fig. 1 for t > 100. 

The DLG strings coarsen with time until well de- 
fined, relatively-narrow longitudinal (i.e., directed along 
x) stripes are formed. (For periodic boundary condi- 
tions, the case of our simulations, each stripe forms a 
ring.) This results into a multi-stripe state, as illustrated 
in figs. 1 and 2. The ordering time in the DLG, de- 
fined as the time the system needs to form the stripes, 
scales with the system size in the direction of the field, 
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FIG. 1: (a) A series of MC snapshots comparing very early 
patterns for the DLG (with an infinite horizontal field) and 
for the standard lattice gas, i.e., zero field (LG) at the same 
time. This corresponds to a 256x256 lattice at T ~ 0.6T<S°. 
The time (in MC steps) is here t = 4 and 10 (from top to 
bottom) for the DLG (left column) and for the LG (right 
column), (b) The same as (a) but at late time, namely, t = 
100, 1,000 and 10,000 (from top to bottom) for the DLG (left 
column) and for the LG (right column). 




FIG. 2: A series of MC snapshots that illustrate (late) growth 
at T ~ Q.GTq' . This is for a rectangular lattice of size L± x 
Ly = 256 x 64 and t = 10 2 , 10 4 , 10 5 , 10 6 , 10 7 and 1,1x10 s 
MC steps, respectively, from left to right. 



L||, since in this case ordered clusters (stripes) percolate 
along the field direction (see below). |26j This is not the 
case for the equilibrium LG, where the ordering time de- 
pends exclusively on system intensive parameters such as 
temperature and density. 

The multi-stripe states are not stable, however. They 
are only partially segregated and, in fact, a definite ten- 
dency towards a fully segregated state with a single stripe 
is generally observed in computer simulations. One may 
also develop simple arguments indicating that, in gen- 
eral, a multi-stripe state will monotonically evolve until 
forming a single stripe. 0,113 It is true that, in prac- 
tice, the complete relaxation may take a very long time. 
More specifically, a macroscopic system may take to de- 
cay into the true stable state a long, macroscopic time 
interval, namely, a time that may show up as mathemat- 
ically infinite in some time scales. In fact, the complete 
relaxation time is observed to increase with system size, 
as first demonstrated in j^. It should also be remarked 
that this property is not a nonequilibrium feature but 
occurs already in the equilibrium (E = 0) case; see, for 
instance, 0, Q and references therein. Slow relaxation 
is a consequence of the conservation of particle density 
p implied by the particle-hole exchange dynamics; this 
induces scale invariance, namely, slow (power-law) evo- 
lution of correlations so that, once enough order sets in, 
all but very small pattern modifications during a single 
MC step are precluded. Consequently, certain individ- 
ual runs sometimes block for a long time in a state with 
several stripes; however, this does not correspond to the 
average behavior. As illustrated by fig. 2, which shows a 
typical evolution, and demonstrated below by our aver- 
ages corresponding to thousand evolutions, the number 
of stripes monotonically decreases with time (see also §}, 
and the whole relaxation can easily be observed in com- 
puter simulations if one waits long enough. 

We next attempt a theoretical description of the relax- 
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ation process. Our interest is on the anisotropic spinodal 
decomposition by which the earliest state with many well- 
defined stripes decays into a single stripe. We shall assume 
that relaxation is a consequence of monomer events caus- 
ing effective diffusion of liquid stripes. (Note that as- 
suming gas stripes here would be completely equivalent.) 
That is, due to single particle processes, liquid stripes 
move transversely as a whole, and may collide and even- 
tually coalesce with one of the neighboring stripes; see 
the late evolution depicted in fig. 2. We notice that co- 
alescence implies evaporation of the gas stripe between 
the two involved liquid stripes. Therefore, given the par- 
ticle/hole symmetry, our assumption is in a sense equiv- 
alent to assuming that growth is due to evaporation of 
stripes; however, the view adopted here allows for a 
more detailed description below. 

In order to evaluate the implications of stripe effective 
diffusion via monomer events, lets assume that stripes 
are well defined, compact and exhibit a (linear) inter- 
face which is rather flat. This is perfectly justified at 
sufficiently low temperature (the case analyzed in de- 
tail here), [2l] and it might hold more generally, in a 
wide region including the center of the miscibility gap 
but excluding the critical region. Under this assump- 
tion, consider a stripe of mean width £ (t) that con- 
sists of M particles whose coordinates along the trans- 
verse (vertical) direction are yj (t) ; j = 1, . . . , M. We 
characterize the stripe position by its center of masses, 
Y cm (t) = M- 1 J2j Vj (t) ■ 

Let us evaluate the mobility coefficient T>£ = 
N mc {(AY cm ) 2 ) which depends on the stripe width £{t). 
Here N mo is the number of monomer events per unit 
time, and ((AY~ cm ) 2 ) is the mean squared displacement 
of the stripe associated to one of the monomer events. 
We think of two possible types of events, each giving a 
different contribution to T>i : 
(A) Evaporation-condensation of particles and holes in 
the stripe surface. Here particles (holes) at the stripe 
interface evaporate to the hole (particle) gas, and con- 
densate later at the same interface. The evolution of 
the evaporated particle (hole) in the bulk can be seen 
as a one dimensional random walk with two absorbing 
walls, the left and right interfaces, respectively. Accord- 
ing to standard random walk theory. |31| the evaporated 
particle (hole) will go again with unit probability to one 
of the (two) possible interfaces, Moreover, the random 
walker will stick again to its original interface with high 
probability, so trapping a particle (hole) from the oppo- 
site interface is unlikely. Consequently, in this case (A), 
N me is simply the evaporation rate. That is, -/V mCj A = 
v exp(— 2T~ 1 Aj), where v is the a priori frequency, 
the sum is over the surface particles and Aj is the num- 
ber of resulting broken bonds. For a flat linear interface, 
particles can only jump transversely away the surface, v 
equals the inverse of the lattice coordination number, q, 
and one may write N me ^A ~ 4g _1 L|| exp(— 2A/T) where 



A is the mean number of broken bonds per evaporation 
event. We multiplied here by 2 to take into account 
evaporation of surface holes that travel within the stripe 
to reach the (same) surface again. On the other hand, 
evaporation processes induce changes AY cm = M~ 1 5y 7 
where Sy is the net particle (transverse) displacement, 
and M « Lm x £(t) for compact stripes. Therefore, 

V[ A) ~ 4<T W) e- 2A/T LjH-*. (2) 

(B) A hole jumps one lattice spacing away within the 
stripe interior. This induces AY cm = 1/M or 0, de- 
pending on the jump direction. One may write -/V m0j B = 
2vph (T) Liilph (T) , where ph is the density of holes, L\\l 
is the volume or total number of sites within the liquid 
stripe, and Ph is the jumping probability per unit time. 
The factor 2 here comes from the fact that a hole mod- 
ifies Y cm when jumping to any of the two directions ±y. 
At low T, ph is small; holes are then rather isolated from 
each other, so that jumps do not modify the number of 
broken bonds, and ph ~ 1. It ensues 

Vf > ~ 2q- l p h L-H- 1 . (3) 

Note that a different dependence of (J2J) and (@ on t is a 
consequence of the fact that the rates N mc ^A and iVme.B 
involve processes consisting of evaporation on the line 
and difussion on the bulk, respectively. 

For (0=5, one has on the average stripes of width £ 
that are separated a distance I from each other. There- 
fore, a given stripe takes a mean time T£ = £ 2 /T>e to 
find (and thus to coalesce with) another one, and this 
causes its width to increase by A£ = I. Consequently, 
d£/ dt ~ A£ t^ 1 = Vii- 1 . Together with © and 
respectively, this implies that mechanism A is character- 
ized by a power law £ ~ i 1 / 4 , and that mechanism B is 
to be associated with I ~ t 1 ^ 3 . Furthermore, assuming 
that pattern growth in the DLG is the result of compe- 
tition between the two mechanisms, and that they are 
independent of each other, T>i = + it follows 

that 

„ .1 ( ^ + u) 

dt L|, V P £ 2 J ' 1 ' 

where cha = 4z/(<5y 2 ) e _2A / T and qb = 2vph. This is our 
general result for the DLG as far as the field E is large, 
e.g., infinite, and the temperature T is low enough so that 
the interfaces, and mechanisms A and B, are sufficiently 
simple as assumed. This is to be compared with the 
Lifshitz-Slyozov- Wagner behavior dt/ At <~ £~ 2 which 
assumes spatial isotropy and diffusion directly governed 
by surface tension. Formally, Q is similar to an equa- 
tion obtained before by assuming isotropic conditions; 
see §1.[H 

The consequences of ar e as follows. Both © 

and imply independently that 

i-wfiv^mf*. (5) 
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The difference is that 9 = and ip = 4 from J5J) while 
one obtains 6 — ocb and ip = 3 from On the other 
hand, for sufficiently late times, I becomes large and 
equation J3J simply solves into 

£(t) ~ at 1/3 + C, (6) 

where a 3 = Sa^L^ 1 and £ = aA/2ae- That is, the pre- 
diction is that hole diffusion within the stripe (mech- 
anism B) will be dominant at late times. A different 
hypothesis, based on the stripe evaporation picture, was 
shown in to imply £ ~ (t/L^ 1 ^ 3 . This coincidence 
is not surprising since, as argued above, the coalescence 
of two particles stripes implies the evaporation of the 
intermediate hole stripe and due to the particle/hole 
simmetry in our system, both mechanisms (stripe diffu- 
sion/coalescence and stripe evaporation) correspond to 
the same physical process yielding the same behavior. 
In order to uncover the close analogy between the two 
pictures, one may notice that, to evaporate a particle 
stripe, many of its particles must cross the surrounding 
hole stripes and stick on the neighboring particle stripes 
(this is so since the particle density in the gas phase re- 
mains almost constant). This particle migration process 
through the surrounding hole stripes is in fact what we 
have called 'hole diffusion within the stripe' in the pres- 
ence of particle/hole simmetry. Hence the fundamental 
mechanism involved in a stripe evaporation is the diffu- 
sion of its constituents through the neighboring stripes. 
This observation is a key one to understand the relation 
between stripe's evaporation and hole (particle) diffusion. 

The effect of mechanism A — surface evaporation and 
subsequent condensation — on growth is more subtle. In 
fact, our theory predicts a crossover from the t 1 / 4 regime 
to the t 1 / 3 regime as time is increased. That is, the two 
mechanisms will have a comparable influence at t ~ r cross 
with 

(4<*a) 3 m 

'cross — 4 -^|| ■ \l ) 

(3ae) 

For times t < t ctoss , mechanism A is dominant and the 
t 1 / 4 behavior is expected, while mechanism B is domi- 
nant for t > T cross and the asymptotic t 1 / 3 growth law 
is then observed. The crossover time r cross is a macro- 
scopic, observable time. Further, we may define the time 
t ss at which a single stripe is reached by the condition 
that £(t) w \L±. One obtains 

L \\ ,h\ (L 2 : 2 

_ 8C 3 [m M2C+|^_ ln 2^ ]} 

L \\ L \\ 

Hence our system is characterized by two different time 
scales, namely, r cross and r ss . They depend on system 
size in a different way. For large systems one generally 



obtains r ss r cross , so that the system converges, after 
a short, perhaps unobservable transient time, to the rel- 
evant i 1 / 3 behavior. However, there are small systems 
for which r ss < r cross . These systems will reach the sta- 
tionary state (a single stripe) before having time to enter 
into the asymptotic t 1 ^ 3 regime. For these small systems, 
the only relevant behavior is the i 1 / 4 one. Therefore, 
there is a size crossover between i 1 / 4 asymptotic behav- 
ior for small systems and t 1 ^ 3 asymptotic behavior for 
large ones. The condition r cross (T, L\\) — t ss (T, L\\, L±) 
defines the crossover size. 

Consider now the parameter 7 = 
Tcross (T,L\\) /t ss (T,L±,L\\) . It follows that the 
i 1 / 3 behavior is dominant for 7 <C 1. However, one 
also has that 7 (T, L±,Lh) — ► for finite T in the 
thermodynamic limit (Lj_,L|| — > 00, L±/L\\ = const.). 
Consequently, the i 1 / 3 growth law is the general one, 
namely, the only one we should expect to observe in a 
macroscopic system. Corrections to this should only 
occur at early times in small systems. This is fully 
confirmed below. 

One may also define a longitudinal length, [2(| [2t| say 
in ~ f a u j where one expects an > 1/3 (given that the 
growth is more rapid longitudinally than transversely). 
This length is only relevant during the initial regime, un- 
til stripes become well-defined, all of them extending the 
whole length Lu . This condition may be taken as defining 
the onset of the multi-stripe state, which may be char- 
acterized by £|| (r ms ) = Lu, from where it follows that 

T ms ~ L^" 11 . Interesting enough, this is on the macro- 
scopic time scale, as for both r cross and r ss (more pre- 
cisely, T cross ~ Lu and t ss ~ LnLjJ. The fact that all 
these relevant times arc on the macroscopic, observable 
time scale confirms that, as argued above, the single- 
stripe (and not the multi-stripe) state is the only stable 
one in general. It is also to be remarked that, once the 
multi-stripe state sets in, the only relevant length is the 
transverse one, £. Of course, this is compatible with the 
possible existence of two correlation lengths describing 
thermal fluctuations at criticality. 

In order to test our predictions, several measures of the 
relevant length in computer simulations were monitored, 
namely: 

• the maximum width of the stripe, £ max , averaged 
over all stripes in the configuration. This maximun 
width is defined as the distance in the direction 
perpendicular to the field between the leftmost and 
the rightmost particles within the stripe; 

• £m = M/L\\, where M = M (t) is the mass, or 
number of particles belonging to the stripe, aver- 
aged over all stripes in the configuration. This mass 
width is defined as the width of a perfectly dense 
stripe with M particles. 



FIG. 3: Time evolution of the relevant length, £(t), as ob- 
tained by different methods, namely, from the number N s of 
stripes (dashed line), from the maximum width, i?max (□) 
from the mass, £m (A), and from the peak of the structure 
function, Is = 27r/fcx,max (O); these quantities are defined in 
the main text. The graphs here correspond to an average over 
600 independent runs for the 128x128 lattice. 



• £ s = L±/2N S , where N s is the number of stripes in 
the configuration. 

After averaging over many independent evolutions, all 
these quantities happen to behave similarly with time. 
Further measures of the relevant length that we define in 
the next section behave in the same way. We shall refer 
to this common behavior, which is illustrated in fig. 3, as 
£ it) . (It is noticeable that, before showing a common be- 
havior, Fig. 3 reveals some significant differences between 
our measures of £ (t) at early times. This confirms the 
more difficult description — not attempted here — which 
is required by the initial regime.) 

In fig. 4 we illustrate our analysis and main results con- 
cerning the (late) time evolution of I it) . The predictions 
above are confirmed and, in particular, "small" lattices 
— fig. 4(a) — happen to behave differently than "large" 
lattices — fig. 4(b). In both cases we plotted £ (t) ver- 
sus t a± for varying aj_, looking for the best linear fit 
£(t) — at a± + C, excluding the initial time regime. The 
upper insets in the figures show the chi square function 
associated to each fit, namely, 



x 2 K) 



v 

E 

i=i 



[£(U)-(at^+C)Y 



at - - 



c 



(9) 



for a least-squares fit to r\ data points using parameters 
a±, a and £. The main graphs confirm the existence 
of a common behavior for all the monitored measures 
of £(t) (indicated by different symbols). These graphs 
also demonstrate that £{t) = at a± + (, with small (, 
during the whole time regime of consideration. On the 
other hand, the upper insets indicate that a±_ is very close 
to j for "small" systems (in fact, for L± < 128) while 




FIG. 4: (a) The main graph shows £ (t) = £ s (t) versus t a± for 
ax = 1/4 in the case of the "small" 64x64 lattice. A similar 
behavior is obtained for any of the studied measures of £ (see 
the main text for definitions), which are represented in the 
insets by different symbols, namely, ^max (empty squares), £s 
(empty circles), and £m (triangles). The upper inset shows the 
chi square function for varying a± as obtained from a series of 
fits; a well-defined minimum is exhibited indicating that a± ~ 
1/4 in this case. The lower inset shows the effective exponent, 
dlog 2 i? / dlog 2 t, as a function of l/£(t); this extrapolates 
to the same value of a±. (b) Same as (a) but demonstrating 
that a± = 1/3 for the "large" L± x L« — 256a;64 lattice (one 
obtains a similar result for larger L||).medskip 



a± ~ | as the system becomes larger, say L± > 256 that 
corresponds to a "large" lattice according to familiar MC 
standards. As an alternative method to analyze I (t) , one 
may evaluate 



o(t) 



diog„ e(t) 

dlog„t 



(10) 



Our prediction is that a(t) = a± — (a±/£(t) , i.e., this 
should provide the exponent a± by extrapolating to large 
£(t) (late time). The insets at the bottom of each fig. 4 
show the results for n = 2. They are in agreement with 
the other method, and again confirm our predictions. 

As indicated above, the size crossover between the i 1 / 4 
and i 1 / 3 asymptotic regimes is expected for a system 
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FIG. 5: The parameter 7 = rcross (T, Ly )/r ss (T, L ± , L\\ ), 
with the characteristic times T cr oss and r ss defined in the main 
text, as a function of L± for L11 = 64. This confirms our dis- 
tinction between "small" and "large" lattices, as explained in 
the main text. 



size such that r cross (T, Ly) = r ss (T, Ly) . 

In order to make this condition explicit, we need to esti- 
mate the amplitudes a a and a-Q m Q; see equations 
Q) and ©. These amplitudes, which state the rela- 
tive importance of surface evaporation/condensation ver- 
sus bulk hole-diffusion, are given respectively by q;a = 
4g~ 1 (<5y 2 ) e~ 2A / T and qb = ^q^ 1 ph- We note that, for 
a sufficiently flat interface (i.e., one that involves micro- 
scopic -but not macroscopic- roughness), (Sy 2 ) ~ O (I) 
and A ~ 5. On the other hand, the excess energy associ- 
ated with an isolated hole is 16, so that pu ~ exp (— I6/T) 
is a rough estimate of the hole density. As depicted in 
fig. 5, it follows numerically, in full agreement with our 
observations, that a± — i is to be observed only at 
early times, earlier for larger systems; to be more spe- 
cific, the crossover for Ly = 64, for instance, is predicted 
for L± ~ 140, which confirms the above; see also figures 
4. 

This behavior may be understood on simple grounds. 
The surface/volume ratio is large initially and, conse- 
quently, mechanism A (based on surface events) is then 
dominant. This is more dramatic the smaller the system 
is. That is, the surface is negligible for macroscopic sys- 
tems, in general, and, as illustrated in fig. 2, even if the 
surface is relevant at very early times, its ratio to the vol- 
ume will monotonically decrease with time. This causes 
hole diffusion in the bulk (mechanism B) to become dom- 
inant, more rapidly for larger systems, as the liquid phase 
is tr ying to exhibit only two surfaces. On the other hand, 
ref. |26J studies the stripe coarsening process in the in- 
finitely driven lattice gas. Pure t 1 / 3 behavior is reported 
assuming the stripe evaporation mechanism. This result 
is perfectly compatible with our results, given that the 



systems in ref. |26| correspond to very large values of L± 
(800 and 960) and small values of Ly (8, 16 and 32). For 
these shapes our theory also predicts the (simple) t 1 / 3 
asymptotic behavior. 



CORRELATIONS AND THE STRUCTURE 
FACTOR 

Consider now the Fourier transform of the pair cor- 
relation function C{x,y]t) = (no.o (t) n x , y (t)), where 



stands for the occupation variable at lattice site 
r— (x,y) . This is the so-called structure factor, S(k,t), 
where k = (fcy, k±). Given that the k\\ dependence is only 
relevant at early times, before the multi-stripe state sets 
in, i.e., for t < r ms , we shall set fc|| = in the following. 
That is, our interest here is on 



S(k X ;t) 



1 



L\\L ; 



^nx, y (t) exp [ikxy] 



(11) 



As illustrated in fig. 6, this function develops a peak at 
kj_ = fcmax (t) immediately after quenching. The peak 
then monotonically shifts towards smaller wave numbers 
with increasing t; in fact, one expects k± — > as t — > 00 
in a macroscopic system. The wave length is = 27r/fc max 
turns out to be an excellent characterization of the rel- 
evant order, namely, it measures both the stripe width 
and the stripe separation during phase segregation. In 
particular, we confirm that is (t) has the common be- 
havior discussed above for length i (t) ; see figs. 3 and 
4. 

The fact that the DLG shows a unique time- dependent 
relevant length, £± = i (t) , has some important conse- 
quences. For example, extrapolating from the equilib- 
rium case (see §f),jjj| one should probably expect dy- 
namical scaling, i.e. 



S(k x ;t) <x £ (t) F [k ± i (t)} 



(12) 



for the anisotropic DLG in two dimensions. This is in- 
deed observed to hold during most of the relaxation and, 
in particular, during all the segregation process after 
formation of well-defined stripes. This is illustrated in 
fig. 7 depicting the scaling function F. A time-dependent 
mean-field model of a binary mixture in shear flow has re- 
cently been demonstrated to exhibit a similar property, 
though involvin g tw o lengths both behaving differently 
from i(t) above. [13 

The structure factor may be obtained by scattering, 
which makes it an important tool in many studies. Ana- 
lyzing further the details of functions S (k±;t) and F (x) 
or, alternatively, the universal function $ (x) = S/£L, 
as observed in computer simulations is therefore of great 
interest (the extra L factor in the definition of $ (x) is 
our finite size scaling ansatz). Experimental studies of- 
ten refer to the mean 'radius of gyration' of the grains as 
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FIG. 6: Time development of the structure factor S(k±;t), 
as defined in the main text, for a "large" lattice L± x Lu = 
256 x 256 during early and intermediate phase segregation. A 
peak grows with time as it shifts towards the small values of 
k ± . 



FIG. 7: The scaling with both time and size of the structure 
function to show that = S {k±; t) /IL, with k — k±£L~ 1 , 
is well-defined and universal, i.e., the same at any time (ex- 
cluding some early evolution) and for any square lattice of side 
L. This plot includes all data for t > 10 4 MC steps and 64x64, 
128x128, and 256x256 lattices. The broken lines illustrate the 
different kinds of behavior of $ (k) that are discussed in the 
main text. 



the slope of the straight portion in a plot of In [S (k, £)] 
versus k 2 .^S^ We checked the validity under anisotropic 
conditions of this concept, which is in fact quite useful 
in equilibrium even outside the domain of validity of its 
approximations. We confirm that S (k±; t) exhibits the 
Guinier Gaussian peak, namely, 



$ (x) ~ exp 



-const. (x 



(13) 



around the maximum x max . More intriguing is the behav- 
ior of $ (x) before the peak, x < x max . Fig. 7 indicates 
that scaling does not hold in this region even at the end 
of our (otherwise long enough) simulations. This is so 
because $ (x) goes as p 2 L/£(t) at k± — 0, and thus de- 
pends on time for very small values of k, breaking the 
scaling observed for larger values of n. However, a de- 
tailed study of data reveals that the scaling function near 
the origin tends with time towards a common envelope 
$ (x) ~ x 1+1 / 3 for x < x < x max ; we do not have 
a simple explanation of this. In any case, this behavior 
breaks down close to the origin, x < xo, where $ (x) — > 
as x — > and t — > oo for the infinite system. 

The behavior after the peak, x > x max , may be pre- 
dicted on simple grounds. The (sphericalized) struc- 
ture factor for (equilibrium) isotropic binary mixtures 
is known to satisfy the Porod's law, S ~ k~^ d+1 ^ at 
large enough k, where d is the system dimension, i.e., 
S ~ fc -3 in two dimensions. The main contribution to 
the large-A; tail comes from the short-distance behavior of 
C(x,y;t). That is, the Porod's region for the DLG may 
be taken to correspond to Aj_ <C fcj 1 <C I (t) , where Aj_ 
stands for a (transverse) thermal length that character- 



izes the smallest, thermal fluctuations. Let two points, ro 
and ro+r, r — (x, y). For any x such that A_l <§C x <C i{t), 
one roughly has that the product n? (t) n?(t) equals +1 
if the two points are on the stripe, and otherwise, i.e., 
if cither an interface exists between them or else the two 
points belong to the gas between stripes. Since x <gC £(t), 
the probability that r crosses more than one interface is 
negligible. For a half-filled system, the probability that 
r*o lies at a particle stripe is ^, and the probability that 
both r*o and ro + r belong to the same stripe is roughly 
\ {£ (t) ~ x) /I (t) . Hence, 



x<£(t). 



(14) 



By power counting, this implies the anisotropic Porod 
law (in two dimensions): 



S(k±;t) 



1 



£{t) k 2 ± 



A ± <fcJ 1 <^(i). (15) 



Therefore, $ (x) ~ x _2 L~ 1 , which is confirmed in fig. 7. 
This is in contrast with the (isotropic) Porod's result. 
The difference is a consequence of the fact that the DLG 
clusters are stripes that percolate in the direction of 
the field, instead of the isotropic clusters of the LG. 
The short-distance pair correlation function for the lat- 
ter is C(r; t) ~ h (1 — \r\/£ (t)), from which one has that 



$(x) 



,-3 



It follows that anisotropy may easily be 



detected by looking at the tail of the structure factor. 

The detailed analysis of S (k±;t) also reveals that, as 
Lu is increased in computer simulations, the anisotropic 
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behavior (f> ~ yt~ 2 crosses over to $ ~ >f~ 3 for larger n\ 
see fig. 7. We believe this reflects the existence of stan- 
dard thermal fluctuations. That is, very small clusters 
of particles occur in the gas in the asymptotic regime 
whose typical size in the direction perpendicular to the 
field is of order Aj_. These very-small asymptotic clus- 
ters are rather isotropic, namely, they do not differ essen- 
tially from the corresponding ones in equilibrium binary 
mixtures. More specifically, for x ~ Aj_, one may ap- 
proximate C(r; t) ~ 1 — \r\/X±(t), which implies the x~ 3 
power-law tail for large x. On the other hand, according 
to (JSJ, the mean stripe width grows as i{t) ~ (t/L^) a 
with a — 1/4 or a = 1/3, depending on the value of 
Lj_. Therefore, the number of stripes at time t is pro- 
portional to L±L^/t a and, for a given time, the number 
of stripes increases with Ly as LJf. We also know that, 
at a given time, the number of small, fluctuating clusters 
is proportional to Lu. Hence the relative importance of 
small clusters due to thermal fluctuations as compared 
to stripes is proportional to Lm. In fact, the xr -3 tail is 
observed for large enough values of L\\ but not for small 
lattices. 



A CONTINUUM DESCRIPTION 

The rigorous derivation of a general continuum analog 
of the driven lattice gas is an open problem. |2JJ Recent 
studies led to the following proposal for a coarse-grained 
density, (f>(r,t) : |34| 

d t <j>(v,t) =T ± Vi0-Vi0+^Vi0 3 +T||Vf0 + V ± £(r,i). 

6 11 

(16) 

Here, the last term stands for a conserved Gaussian 
noise representing the fast degrees of freedom, and t±, 
tii and A are model parameters. Compared to previous 
proposals, [3113 this Langevin type of equation amounts 
to neglect a nonlinear current term, — aV\\(f) 2 , that was 
believed to be essential (relevant) at criticality. How- 
ever, one may show that, at least in the limit E — > oo, 
the coefficient a cancels out (due in this case to a sub- 
tle saturation effect). [34| In fact, recent scaling analysis 
has unambiguously confirmed that a particle current is 
not relevant and that equation (ll(jl) captures the correct 
critical behavior of the DLG.|27J, |22j Consequently, an 
important question is now whether (|16|) reproduces also 
the kinetic behavior of the DLG as described in previous 
sections. We present in this paper a first confirmation 
that, as compared with other approaches, 24] (|16(l is in- 
deed a proper continuum description of the DLG kinetic 
relaxation. A more complete study of the kinetic con- 
sequences of this equation will be presented elsewhere. 

In order to numerically integrate l|16[l , let us introduce 
the indexes i,j = l,...,Nto represent, respectively, the 




FIG. 8: Series of snapshots as obtained from equation 1171 
for the 256 x 256 lattice with parameters as given in the main 
text. Time (arbitrary units) is t = 0, 10, 100, 200, 500, and 
1000 respectively, from left to right and from top to bottom. 



two components of r = (x±, Xn) . One thus makes a triv- 
ial discretization of th e sp ace, and then of the time by 
Cauchy-Euler method. |36j The result is a set of N 2 — 1 
coupled nonlinear equations, namely, 

<b{i,j;t + At) = <Kz,j,i) + Ai[T ± Vi0-Vl>+^Vi0 3 

+T 1 |Vjjfl + VAtV. L £(*,j;<)- (17) 

This equation is to be solved by the com- 
puter. With this aim, we may write Vj_£(i,j;£) = 
[£(i+l,j;t)-£(i-l,j;*)]/2Ax ± and </>{i,j;t) = 
(j>(iAxx,jAx\\;t) where Ax± = L±/N and Axy =L\\/N. 
The maximum value of Ax± is thus limited by the in- 
terface width. For fig. 8, which concerns a 256 x 256 
lattice (N — 256) we — rather arbitrarily — used 
Ax± = Axn = 1.7, and At — 0.05, which produce a 
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locally stable solution. The parameters t±, tu and A, 
are fixed on the basis of its physical meaning. The mass 
terms n and t± represent temperatures along the lon- 
gitudinal and transverse directions, respectively, relative 
to the critical temperature, i.e., t± ~ (Tj_ — Tg?). Given 
the anisotropy of phase segregation, with longitudinal 
interfaces only, r± < and tii > 0. On the other hand, 
\t±\ should be small enough to allow for a relatively fast 
evolution. Our choices for fig. 8 are t± = —0.25, t\\ = 0.5 
and A = 1. 

It is remarkable that, in spite of some apparent sim- 
ilarity, the problem here differs from the one in the 
study of (standard) spinodal decomposition by means 
of the isotropic (E = 0) Cahn-Hilliard equation. In 
equilibrium, [2Z| one usually assumes that the influence 
of noise on growth, which is then assumed to be directly 
driven by surface tension, is negligible far from criticality. 
The noise term in (|17ll may be expected to be important 
in a more general context, however. That is, as described 
in Sfc the DLG develops striped patterns in which surface 
tension smooths the interfaces but has no other dominant 
role on the basic kinetic events. Consequently, neglect- 
ing the noise in l|17fl would turn metastable any striped 
geometry after coarsening of strings, which is not accept- 
able (see 

Finally, it is interesting to notice that if a one- 
dimensional structure is assumed, and the gradient in 
the direction parallel to the field in Eq. [TJ^is eliminated, 
then this equation reduces to the one-dimensional time- 
dependent Ginzburg-Landau model in Q . There it was 
found a ln(t) growth at zero temperature and a crossover 
from ln(t) to t 1 / 3, at finite temperatures. 

CONCLUSION 

This paper presents a theoretical description of spin- 
odal decomposition in the DLG, and compares it with 
new data from a kinetic Monte Carlo study. This is also 
compared with the kinetic implications of a Langevin, 
continuum equation that had previously been shown to 
capture correctly the critical behavior of the DLG. The 
resulting picture from these three approaches, which is 
summarized below, should probably hold for a class of 
higly-anisotropic phase segregation phenomena. In fact, 
our results provide a method for analyzing experiments 
that could be checked against laboratory realizations of 
the DLG, i.e., the case of phase segregation under biased 
fields or other influences such as electric fields, gravity 
and elastic stresses. 

Immediately after a deep quench, there is an early 
regime in which anisotropic grains develop. They tend to 
coarsen to form small strings that then combine into well- 
defined thin stripes. Such nucleation and early coarsen- 
ing (fig.l) seem governed by surface tension at the string 
ends competing with other both surface and bulk pro- 



cesses. This complicated situation typically extends less 
than 10 3 MC steps in computer simulations, which cor- 
responds to a very short macroscopic time, so that it 
would be hardly observable in experiments. As a mat- 
ter of fact, most of the system relaxation proceeds by 
coarsening of stripes until full segregation (fig. 2). Sur- 
prisingly enough, this re gime , which has been studied for 
more than a decade now!|22^-[2fij| happens to be theoreti- 
cally simpler than the corresponding one for the isotropic 
case. H|- [Hi 

The evolution from many stripes to a single one mainly 
proceeds by competition of two mechanisms: (A) evap- 
oration of a particle from one stripe surface and sub- 
sequent deposition at the same surface, and (B) diffu- 
sion of a hole within the bulk of the stripe. The first 
one dominates initially (and lasts more for smaller sys- 
tems), when the surface/volume ratio is relatively large. 
Mechanism A implies that the relevant length (as defined 
in fig. 3) increases with time according to lit) ~ i 1 / 4 . 
The surf ace/ volume ratio decreases with time, however, 
and mechanism B soon becomes dominant. This im- 
plies £ (t) ~ t 1 / 3 which is the general prediction for a 
macroscopic system (cf. figures 4 and 5).|2a| This was 
obtained before by assuming coarsening of two (liquid) 
stripe s by evaporation of the gas stripe placed between 
them; (U see also [z| Note that the t 1 ^ 3 law is precisely 
the behavior which is acknowledged to be dominant un- 
der isotropy, but this has a different origin in the equilib- 
rium case. [1,0] Note also that surface tension determines 
evaporation rates but has no other influence on mecha- 
nisms A and B. 

The i 1 / 3 growth law, JJjJl, is perfectly confirmed by 
the DLG data (fig. 4b). This indicates time-scale in- 
variance. In fact, such invariance was demonstrated for 
the isotropic case, in which the situation is somewhat 
more involved (§1). The invariance property may be 
better analyzed by looking at the structure factor trans- 
versely to the drive, S(k±,t) (fig. 6). This exhibits dy- 
namic scaling, i.e., it remains self-similar during phase 
segregation. [39| More specifically, = S (k±;t) /£L, 

with yc = k±£L~ 1 , is universal, namely, the same at any 
(sufficiently late) time t and for any square lattice of side 
L. Furthermore, the function $(x) has a well-defined 
shape. In particular, it exhibits the Guinier Gaussian 
peak, and this is followed by the anisotropic Porod de- 
cay, $ (x) ~ x~ 2 and then by a thermal tail $ (x) ~ x~ 3 
(fig. 7). Also noticeable is the fact that the the param- 
eter to scale along the S axis is J(t) = £ and not £ 2 as 
under isotropy. 

Our results in this paper have two main restrictions, 
both due to the great computational effort required by 
this problem. [2^ Firstly, they follow from an extensive 
analysis of only one phase-diagram point, namely, p = 
E = 00, and T = 0.8TjS. However, our own observa- 
tions (including brief investigation of other points), to- 
gether with an extrapolation of the many results known 
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for the isotropic case, strongly suggest that the picture in 
this paper holds within a larg e domain around the cen- 
ter of the miscibility gap.|40| In fact, the scaled struc- 
ture factor for isotropic systems was shown to be al- 
most independent of density and temperature, and even 
the substance investigated, in a wide region below the 
coexistence line. [To| Our consideration of only a two- 
dimensional system does not seem a real restriction nei- 
ther. That is, adding an extra (transverse) dimension 
should not essentially modify the picture here.|4ll| 

It would be interesting to look next in the labora- 
tory for both time-scale invariance and t 1 ^ 3 growth under 
highly anisotropic conditions. In fact, there are some ev- 
idences of such behavior in sheared fluids (§1), and one 
may think of some more direct experimental realizations 
of the driven lattice gas. In particular, coarsening striped 
patterns very similar to those observed in our system 
are found in some intringuing experiments on granular 
binary mixtures under shaking. |42j We think that the 
mechanisms we propose in this paper should help the 
understanding of such experimental results. In general, 
we hope our observations will motivate both experiments 
and future more complete theories. 
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